OPEN Q ACCESS Freely available online 



•0-PLOS I ONE 



Quantitative Trait Loci for Resistance to the Congenital (S\ 
Nephropathy in Tensin 2-Deficient Mice c2L 

Hayato Sasaki 1 , Nobuya Sasaki 2 *, Tomohiro Nishino 1 , Ken-ichi Nagasaki 3 , Hiroshi Kitamura 4 , 
Daisuke Torigoe 1 , Takashi Agui 1 

1 Laboratory of Laboratory Animal Science and Medicine, Department of Disease Control, Graduate School of Veterinary Medicine, Hokkaido University, Sapporo, Japan, 

2 Laboratory of Laboratory Animal Science and Medicine, Faculty of Veterinary Medicine, Kitasato University, Towada, Japan, 3 Section of Biological Safety Research, 
Chitose Laboratory, Japan Food Research Laboratories, Chitose, Japan, 4 Department of Veterinary Physiology, School of Veterinary Medicine, Rakuno Gakuen University, 
Ebetsu, Japan 



Abstract 

The ICR-derived glomerulonephritis (ICGN) mouse is a chronic kidney disease (CKD) model that is characterized 
histologically by glomerulosclerosis, vascular sclerosis and tubulointerstitial fibrosis, and clinically by proteinuria and 
anemia, which are common symptoms and pathological changes associated with a variety of kidney diseases. Previously, we 
performed a quantitative trait locus (QTL) analysis to identify the causative genes for proteinuria in ICGN mice, and found a 
deletion mutation of the tensin 2 gene (Tns2" ph , MGI no: 2447990). Interestingly, the congenic strain carrying the Tns2 nph 
mutation on a C57BL/6J (B6) genetic background exhibited milder phenotypes than did ICGN mice, indicating the presence 
of several modifier genes controlling the disease phenotype. In this study, to identify the modifier/resistant loci for CKD 
progression in 7hs2-deficient mice, we performed QTL analysis using backcross progenies from susceptible ICGN and 
resistant B6 mice. We identified a significant locus on chromosome (Chr) 2 (LOD = 5.36; 31 cM) and two suggestive loci on 
Chrs 1 0 (LOD = 2.27; 64 cM) and X (LOD = 2.65; 67 cM) with linkage to the severity of tubulointerstitial injury. One significant 
locus on Chr 13 (LOD = 3.49; approximately 14 cM) and one suggestive locus on Chr 2 (LOD = 2.41; 51 cM) were identified as 
QTLs for the severity of glomerulosclerosis. Suggestive locus in BUN was also detected in the same Chr 2 region (LOD = 2.34; 
51 cM). A locus on Chr 2 (36 cM) was significantly linked with HGB (LOD = 4.47) and HCT (LOD = 3.58). Four novel epistatic 
loci controlling either HGB or tubulointerstitial injury were discovered. Further genetic analysis should lead to identification 
of CKD modifier gene(s), aiding early diagnosis and providing novel approaches to the discovery of drugs for the treatment 
and possible prevention of kidney disease. 
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Introduction 

ICR-derived glomerulonephritis (ICGN) mice are an inbred 
strain that exhibit glomerulosclerosis characterized by marked 
glomerular hypertrophy, mesangial expansion, foot process 
effacement of podocytes and thickening of the glomerular 
basement membrane. The ICGN mouse is also a chronic kidney 
disease (CKD) model that presents common symptoms and 
pathological changes associated with a variety of kidney diseases, 
such as albuminuria, hypoproteinemia, hyperlipidemia, anemia 
and systemic edema and eventual end-stage renal failure [1,2]. In 
ICGN mice, nephrotic syndrome occurs at the suckling stage, 
which is unique in spontaneously nephropathic strains of mice or 
rats [3]. Although glomerular alteration is observed from two 
weeks after birth, renal function is preserved until three to four 
months after birth, when ICGN mice begin to show tubulointer- 
stitial fibrosis [4]. This pathological process in ICGN mice 
conforms well with the clinical presentation of glomerular disease 
in humans [5]. 



Proteinuria is intimately involved in the dysfunction of 
podocytes or slit diaphragm, and glomerulosclerosis (GS) sequen- 
tially starts with podocytopenia. A line of evidence suggests that 
podocytes act as the main component of this barrier as mutations 
in a number of podocyte-specific genes (NPHS1, NPHS2, CD2AP, 
ACTJV4, TRPC6, PLCE1, AITH9) responsible for GS have been 
identified [6,7]. Therefore, podocyte loss is a common determining 
factor for progression toward many types of kidney disease [8,9]. 
Previously, we performed a QTL analysis to isolate the causative 
genes for ICGN phenotypes, and identified a major QTL on 
chromosome (Chr) 15 that is identical to a single recessive locus 
causing albuminuria, and found a deletion mutation of the tensin 2 
gene (Tns2, also known as Tend). This mutation (Tns2 n ^ A , MGI no: 
2447990) leads to the loss of tensin 2 in podocytes, parietal 
epithelial cells and tubular epithelial cells [10]. The tensin (Tns) 
family of multidomain scaffold proteins (including Tnsl, Tns2, 
Tns3 and Tns4) bind to the cytoplasmic tail of fi-integrins, and 
localize to adhesions that anchor stress fibers in cells [11]. 
Therefore, Tns2 is thought to be a component of the actin 
cytoskeletal structures linking actin filaments with focal adhesions 
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and plays a role as an intracellular signal transduction mediator 
through integrin [11]. It has been reported that Tns2 may be 
involved in the expression of nephrin (Nphsl) and synaptopodin 
[Sjnpo) in podocytes [12]. However, the mechanism by which Tns2 
deficiency leads to CKD is unclear. 

In humans, the prevalence rates of CKD vary considerably by 
race, sex and age [13]. Multiple studies, such as genome- wide 
association studies, have identified many genetic factors that 
participate in kidney disease, although it is estimated that a vast 
amount of genetic factors for CKD remain to be identified [14]. 
CKD has a complex etiology resulting from the combination of 
multiple genetic and environmental factors, each with small 
effects. Epistatic interactions among susceptibility modifier loci 
make many of the loci difficult to detect in genome-wide 
association studies [15,16]. Several reports using mouse models 
of nephropathy have demonstrated similar phenomena in that the 
effects of nephropathy-inducing factors were attenuated in other 
genetic backgrounds [17,18,19]. Recently, we have shown that 
congenic strains carrying the Tns2 n> ' 1 ' mutation on a B6 or 129x1/ 
SyJ (129) genetic background exhibited milder phenotypes than do 
ICGN mice [20,21]. In contrast, other researchers have reported 
DBA/2J (D2) congenic mice with Tns2 nph were susceptible to CKD 
at a level similar to ICGN mice [22]. This suggests that B6 and 
1 29 mice have resistant genes to CKD induced by TnsT^ 11 . In the 
present study, we performed a genome-wide linkage analysis in 
(B6xICGN) F, xICGN backcross mice carrying Tm2 npk to identify 
the CKD resistance genes and elucidate the mechanism by which 
Tns2 deficiency leads to CKD. 

Materials and Methods 

Animals 

ICGN mice were purchased from Osaka National Institute of 
Biomedical Innovation. C57BL/6J (B6: Charles River) mice were 
backcrossed to ICGN mice. Tns2 (nphl nph) N2 backcross progenies 
were obtained by mating ICGN male mice to Tns2 (+/nph) Fj 
hybrid female mice heterozygous for Tns2" t ' h . Backcross mice 
homozygous for Tns2" pl1 were born in accordance with Mendelian 
inheritance and at an equal sex ratio. h6-Tns2" 1 ' 1 ' congenic mice 
were generated as described previously [20]. The animal room 
was air-conditioned at 22±4°C, maintained at 40-60% relative 
humidity, and mice were maintained under a 12 hr light-dark 
cycle. A standard laboratory diet, CE-2 (Nihon Clea, Tokyo, 
Japan), and tap water were available ad libitum. In this study, 16- 
week-old mice were selected for all analyses as the ICGN mice 
show severe albuminuria and anemia, which are the end stages of 
a variety of pathological conditions. A humane end point was 
applied when the mice with severe anemia became moribund. 

Ethical Statement 

All research and experimental protocols were adhered to the 
Regulation for the Care and Use of Laboratory Animals, 
Hokkaido University and approved by the President of Hokkaido 
University following to the review of the Institutional Animal Care 
and Use Committee (Approval ID: No. 110226). 

Phenotypic Analysis 

Blood samples were collected from the postcava under 
isoflurane anesthesia. Hemoglobin concentration (HGB), hemat- 
ocrit (HCT) and blood urea nitrogen (BUN) were measured as 
described previously [20]. 

Formalin-fixed and paraffin-embedded kidney blocks were 
sectioned at a thickness of 2 (J.m, and stained with Periodic acid- 
Schiff (PAS). In this study, two types of histopathological scores 
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Figure 1. Hemoglobin concentration and BUN and kidney 
injury in backcross mice. (A, left) Distributions of HGB and BUN in 

(ICGN xB6) F, xlCGN backcross progenies. (A, right) Distributions of HGB 
and BUN in B6-Tns2" ph , ICGN, (ICGNxB6) F„ backcross mice. (B) 
Glomerular index and tubular index distributions in (ICGNxB6) F q xlCGN 
backcross progenies. B6con, B6-Tns2 nph ; BC, backcross mice; NS, no 
significant; asterisk, P<0.001. 
doi:1 0.1 371 /journal.pone.0099602.g001 

were adopted to quantitate the severities of renal damage. A 
glomerular index, in which the severities were based on glomerular 
injury, was calculated from the scores for twenty glomeruli 
randomly selected from each sample. The severity of glomerular 
injury was scored as follows: 0, no abnormality or nearly normal; 
1, expansion of the mesangial matrix; 2, severe expansion of the 
mesangial matrix; 3, hyalinosis or collapse (Fig. SI). First, the 
highest and lowest scores were removed to exclude bias in random 
selection. The sum of the eighteen glomerular scores was then 
taken as the glomerular index. However, if there were glomeruli 
that scored 3 points, the index was the sum of 36 and the number 
of glomeruli that scored 3 points. This calculation method means 
that a kidney with hyalinosis or collapse has a more severe 
glomerular injury than a kidney with only mesangial expansion. 
The tubular index was the sum of two histopathological 



parameters, urinary cast and tubule dilation. The severity of 
tubule dilation was performed using scale 0 to 2:0, no abnormality; 
1, mild; 2, severe. One point was then added to the tubular index if 
a urinary cast was detected (Fig. S2). These analyses were 
performed by two independent observers. There were remarkable 
congruences in both blinded scores [the interclass correlation 
coefficients (ICCs) = 0.96]. 

Genotyping Analysis 

We used 246 (male =124, female =122) Tns2 (nph/nph) N2 
backcross mice for a genome-wide scan. Extraction of genomic 
DNA from tail tips was performed by a standard method. Tns2 
genotypes were determined by PCR using genomic DNA from the 
tail tips. TnsT ph homozygous mice were identified using the 
primers AGACACCACCAGCACCTTCT and CGATC- 
CAGCTCCTGTCTTTC, which distinguish the TnsT ph allele 
by the 8-base deletion in exon 18 of Tns2 [10]. For QTL analysis, 
a total of 92 informative microsatellite markers and coat color 
were used for the genotyping analysis, as listed in Table SI. The 
map positions of the microsatellite markers were based on 
information from the Mouse Genome Informatics of Jackson 
Laboratory (MGI; http://www.informatics.jax.org/, MGI_4.41). 
PCR was carried out with a cycling sequence of 95 °C for 5 min 
(one cycle), followed by 35 cycles consisting of denaturation at 
95°C for 30 sec, primer annealing at 55°C for 30 sec, and 
extension at 72°C for 30 sec. Amplified samples were electropho- 
resed with 12% polyacrylamide gels and stained with ethidium 
bromide. The stained gels were then visualized and photographed 
under an ultraviolet lamp. 

QTL Analysis and Statistics 

QTL analysis was performed using the Map Manager QTXb20 
software program that uses a maximum likelihood algorithm with 
"interval mapping" and "simultaneous search", and permits better 
localization of loci and exclusion mapping [23]. Recombination 
frequencies (%) were converted into genetic distance (cM) using 
the Kosambi map function. This program provides linkage data as 
a likelihood ratio statistic (LRS) score. For each Chr, the LRS 
values were calculated by 10000 random permutations of the trait 
values relative to the genotypes of the marker loci. Genome-wide 
significance thresholds were calculated in terms of LRS by 10000 
random permutation tests based on the established guidelines. The 
logarithm of odds (LOD) score was calculated by dividing the LRS 
value by 4.61. LOD values of the suggestive and significant 
thresholds for each trait are summarized in Table 1. Significant 
linkage was determined according to the traditional threshold 
(LOD score >3.3) [24]. Statistical analyses were performed using 
the Stat- View program (SAS Institute Inc. Cary, USA). Spear- 
man's rank correlation analysis was used to find the correlation 
between two sets of tested traits. 

Results 

7. Phenotypic Characterization of Tns2 nph Backcross Mice 

To verify the effect of genetic background on susceptibility to 
CKD induced by TnsT ph , we produced a total of 246 (ICGNxB6) 
Fi xlCGN backcross mice and performed blood testing and 
histopathological analysis at the beginning of a genome-wide 
linkage analysis. Blood samples collected at 16 weeks old were 
analyzed. In CKD, a high BUN level and low HGB signify renal 
impairment and renal anemia, respectively, which is also a 
common symptom in ICGN mice (Fig. 1A, right). In backcross 
progenies, wide spectrums of BUN and HGB were observed from 
the various degrees of B6 type to that of ICGN types (Fig. 1A). 
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Figure 2. Linkage of nephropathy to chromosome 2. (A) Genome-wide linkage analysis of the tubular index. (B) LOD plots show the linkage of 
nephropathic traits to chromosome 2. HGB (solid line) and tubular index (dotted line) yielded significant LOD scores. The approximate 95% 
confidential intervals: 19.6-56.4 cM (tubule), 26-73 cM (glomeruli), 26-55 cM (HGB), 23-59 cM (HCT), 27-72 cM (BUN). (C) LOD curve and bootstrap 
histogram of the QTL on chromosome 2 for tubular index. 
doi:1 0.1 371 /journal.pone.0099602.g002 



These results imply the existence of multiple modifier loci for 
GKD caused by Tns2 nph . In addition, the distribution of BUN in 
backcross males, as compared with females, was statistically more 
severe (P<0.001, Mann-Whitney's U test). 

Figure IB shows the box plot for glomerular index and the 
histogram for tubular index. The scores histologically evaluate 



glomerular injury and tubulointerstitial injury, respectively. The 
backcross mice showed wide range of the histopathological 
severity, from B6 types to ICGN types (see also Fig. SI, S2, S3). 
There were also significant sexual dimorphism in both scores, with 
males exhibiting a higher level of severity than females (glomerular 
index P<0.001, tubular index P<0.05, Mann- Whitney's U test). 
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Figure 3. LOD score curves on chromosome 13. (A) A significant QTL for glomerular index was detected on chromosome 13. The approximate 
95% confidential intervals: 0-31.5 cM. (B) Glomerular index at the peak LOD score (D13Mit60). The resistance phenotype associated with ICGN/B6 
genotype (P<0.001, Mann-Whitney's U test). 
doi:1 0.1 371 /journal.pone.0099602.g003 



PLOS ONE | www.plosone.org 



4 



June 2014 | Volume 9 | Issue 6 | e99602 



QTL for Resistance to Congenital Nephropathy in Tns2-Deficient Mice 



a 1 

3 



- ICGN r ICGhl 

- ICGN J BG 




ICGN / ICGN ICGN / B6 

O2Mk380 



14 



?13 J 

CD 

S.J 



11 




D19MA33 

- ICGN f ICGN 

- ICGN j B6 



ICGN / ICGN 



ICGN t B6 



I 

3 

,3 





r D2MMQ2 

-O- ICGN t ICGN 

V -■- ICGN f B6 


**- 





ICGN / ICGM ICGNB6 

D2MH293 



3 
J3 



D3M>t1?9 

- ICGN i 1 ICGN 

- ICGN 1 66 




D2MH380 



ICGN J ICGN ICGN / BG 

DlOMtf2T1 



Figure 4. Epistasis associated with QTLs. Epistatic interactions were detected associated with the markers on chromosome 2 and 10. Open 
circles and closed circles represent homozygous and heterozygous variants for the tested markers, respectively. Asterisks signify significant 
differences between the two genotypes (**P<0.001; *P<0.01; Mann-Whitney's U test). 
doi:1 0.1 371 /journal.pone.0099602.g004 



These differences were presumably due to the protective effect of 
estrogen in the female glomeruli [25]. 

The symptom severity of all tested traits exhibited significant 
correlations with each other according to Spearman's rank 
correlation analysis (P<0.01). However, the Spearman's rank 
correlation coefficients were comparatively smaller for BUN 
combinations ( | r | = 0.27 ~ 0.48) than for the other sets 
(|r| = 0.4 ~ 0.65). 

2. Genome-wide Linkage Analysis for the Mapping of 
CKD-Resistant Genes in Backcross Mice 

To identify the modifier genes for CKD, we carried out a 
genome-wide linkage analysis using Tns2 n ^ 1 ' backcross mice. HGB, 
HCT, BUN, glomerular index and tubular index at 16 weeks of 
age were used as indicators for CKD. Table 1 summarizes the 
microsatellite markers linked to the phenotype, LOD score, 
interval and phenotypic values of each genotype. The QTL 
analysis of the tubular index, which shows the intensity of 
tubulointerstitial injury, identified one significant locus on Chr 2 
and two suggestive loci on Chrs 10 and X, respectively (Fig. 2A, 
Table 1). The significant QTL had a peak LOD score of 5.36 and 
explained 10% of the variance. Regarding the other nephropathic 
traits (HGB, HCT, BUN and glomerular index), significant or 
suggestive loci were also detected on Chr 2 (Fig. 2B, Table 1). 



Within these QTLs, each trait exhibited CKD resistance 
associated with the ICGN/B6 genotype. 

The QTL analysis of the glomerular index, which shows the 
severity of GS, identified one significant locus (peak LOD = 3.49) 
on Chr 13 (Fig. 3A). Grouping backcross mice by their genotype at 
the marker near the linkage peak, D13MU60 (14.44 cM), ICGN/ 
B6 heterozygous mice showed statistically milder glomerular 
injury than did ICGN/ICGN homozygous mice (P<0.001, 
Mann- Whitney's U test, Fig. 3B). Although the linkage of the 
tubular index to Chr 13 was similar to that of glomerular index, 
the peak LOD score was only 1.73. No suggestive linkages were 
observed for hematological phenotypes on Chr 13 (Fig. 3A). 

Epistatic interactions were searched between all of the tested 
markers by MapManager QTXb20 (P<0.000001). We detected 
three significant interactions in the tubular index associated with 
markers on Chrs 2 and 10, which showed significant or suggestive 
linkages to tubular index (P<0.01, ANOVA test, Fig. 4). Genotype 
variations in the tubular index at D2MU380 (40.89 cM), D2MU293 
(17.24 cM) and D10MU271 (72.31 cM) were only expressed in the 
presence of the ICGN/ICGN genotype at D19MU33 (51.76 cM), 
the ICGN/ICGN genotype at D2MU102 (57.65 cM) and the 
ICGN/B6 genotype sxD3Mitl29 (80.49 cM), respectively. The B6 
allele at D19MU33 also suppressed the genotypic effect of 
D2MU380 in HGB. D2MU293 (17.24 cM), which lies outside the 
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Table 2. Podocyte 


related genes on Chr 2 QTL. 




Candidate genes 


Interacting proteins (XPodNet) 


Enkur 


Pakl, Pik3r2 


Apbblip 


Enah, Vasp 


Pax8 


Paxipl, Wt1 (Chr 2:55 cM) 


Grin 1 


Ctnnbl, Src (Chr 2:78 cM), Actb, Spnb2, Prkce, Camk2a, Camk2b, Dnml (Chr 2:22 cM), Ppp3ca 


Notch! 


Wwpl, Psenl, Psen2, Ep300, Rbpj, Wdr12, Fbxw7, Kat2a, Kat2b, Smarcd3 


Vav2 


Lcp2, Racl, Nckipsd, Ptpn6, Src (Chr 2:78 cM) 


Rxra 


Trim24, Rarb, Rara, Ppara, Thra, Thrb, Zfp423, Nsdl, Med25 


Olfml 


N.D. 


Tscl 


Ezr, Ywhab (Chr 2:85 cM) 


Abll 


Nckl, Pxn 


Lamc3 


Nidi, Nid2, LamaS (Chr 2:103 cM) 


Dnml 


Atp9a (Chr 2:89 cM), Grinl (Chr 2:17 cM), Grin2b, Sh3gl2, Itsnl, Pfn2, Prnp (Chr 2:64 cM), Ctsl, Ap2a2, Kif2a, Kif5a, Tubala 


Eng 


Tgfbrl, Tgfbr2 


Cdk9 


Cdknla, Cdknlb 


Angptl2 


N.D. 


Lmxlb 


Ldbl, Tcf3, Bmpr2 


Mapkap 1 


Aktl 


Crb2 


N.D. 


LrpW 


Pickl 


Zeb2 


Smad2 


Prpf40a 


Wasl 


Pkp4 


Ctnnbl 


Pla2r1 


N.D. 


Itga6 


Lama5 (Chr.2:103 cM), Itgbl, Itgb4, Nisch 


Wipfl 


Cttn 


Hoxdl 


Ldbl (Chr.1 9:38.75 cM) 


Itgav 


Itgbl, Itgb3, Itgb5, Sppl, Plaur, Col4a3 


Ctnndl 


Yes, Cdh2, Cdh4 (Chr.2:102 cM), Cdh22 (Chr.2:85 cM) 


Lrp4 


Pickl 


Syt13 


Nrxnl 


Wtl 


Pawr, Wtip, Fhl2, Fhl3, Pax2 (Chr.1 9:38 cM), Pax8 (Chr.2:16 cM) 


Pax6 


Paxipl 


N.D.: not detected. 

doi:1 0.1 371 /journal.pone.0099602.t002 



95% confidence interval for the QTL of the tubular index (19.6— 
56.4 cM), showed a significant epistatic interaction with 
D2MU102, whereas the inner marker D2MU369 (24.51 cM) did 
not, although both markers exhibited significant LOD scores of 
3.6 and 4.4, respectively. On the other hand, D2MU380 
(40.89 cM) showed a significant epistatic interaction with 
D19MU33, as did D2MU66 (49.45 cM) and D2MU102 (57.65 cM) 
(data not shown). These results suggest multiple QTLs for tubular 
index on Chr 2. The bootstrap analysis also supports this 
hypothesis (Fig. 2C). Next, to identify the candidate genes, we 
screened the QTLs on Chrs 2 and 1 3 for podocyte-related genes 
using XPodNet, and protein-protein interactions in the podocytes 
[19]. In total, 46 candidate genes are listed (Table 2 and 3). 

We also detected seven genetic loci associated with blood cell 
traits (Table S2). These loci (except for Chr 2) were not linked to 
the histopathological traits, so that the loci only influenced 
hematological parameters between the two mice, independent of 
CKD. Correspondingly, numerous genetic loci of this kind have 



been identified in recent genome-wide linkage analyses of humans 
and other mouse strains [26,27,28,29]. 

Discussion 

In this study, we report the mapping of genetic loci responsible 
for the resistance to congenital nephropathy caused by Tns2 
deficiency in a B6 genetic background. Analysis of backcrossing 
between resistant B6 mice and susceptible ICGN mice identified 
significant QTLs on Chrs 2 and 13. Among the tested traits, 
tubulointerstitial injury (tubular index) and anemia (HGB and 
HCT) were strongly associated with the same locus on Chr 2. In 
patients with CKD, anemia develops due to decreased renal 
erythropoietin synthesis resulting from tubulointerstitial damage. 
A significant locus for tubular index, ranging approximately from 
19.6 to 56.4 cM on Chr 2, which overlaps the QTLs of the other 
tested traits, appears to be separated into three QTLs with peaks at 
22, 31 and 50 cM (Fig. 2). The proximal QTL (22 cM) has 
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Table 3. Podocyte-related genes 


on Chr 13 QTL. 




Candidate genes 


Interacting proteins (XPodNet) 


Net) 


Magil 


Actn2 


Camk2a 


Nid) 


Lama2, Lamcl, Lamc2, Lamc3, Hspg2, Eln 


E2f3 


Ccndl 


Agtrla 


Agtrap 


Ripk) 


Atg5 


Nml 


N.D. 


Bmp6 


Bmprla 


Nedd9 


Src, Nck2, Ptk2 


Gadd45g 


Cdknla 


Nsdl 


Rxra 


Slc34a1 


Itgav 


Dbn) 


Gjal, Lrrk2 


Ntrk2 


Cdk5, Cdk5r1 


N.D.: not detected. 

doi:1 0.1 371 /journal.pone.0099602.t003 



epistatic interactions with D2MU102 and D2MU282 (Fig. 4, data 
not shown). The region between the two markers exhibiting a peak 
LOD score of 3.8 in tubular index includes the syntenic locus for 
albuminuria previously reported in humans, rats and mice [30,3 1] . 
The middle QTL (3 1 cM) maps to the same region that was found 
in a (B6 xD2) F 2 intercross and in a (B6 xA/J) F 2 intercross as a 
second QTL for albuminuria on Chr 2 [30,31]. The distal QTL 
(50 cM), where the glomerular score and BUN showed a 
suggestive linkage, maps to the same region as the suggestive 
locus for proteinuria that was found in an HIV-associated 
nephropathy (HIVAN) model mouse on a mixed genetic 
background of B6 and FVB [32]. The syntenic human region 
(2q32. 1) is associated with the early onset of end-stage renal disease 
in black families enriched for nondiabetic nephropathy [33]. 
Taken together, the resistance to CKD in Tns2-deficient B6 mice 
can be mainly exerted in combination with the multiple loci on 
Chr 2. 

While Tns2 mRNA is broadly expressed in the kidney and other 
organs, significant abnormalities (i.e., podocytopenia) were only 
observed in the kidney of ICGN mice [1]. In the glomeruli, Tns2 
expression is localized in podocytes [10,34]. Therefore, the 
dysfunction of Tns2 in podocytes has been speculated to cause 
renal dysfunction. Tns is an actin-binding protein that is present in 
focal adhesions and has been suggested to play an important role 
in the actin cytoskeleton and integrin-mediated signaling [11]. 
Recently, it has been reported that Tns2 may bind to integrin (3 1 , 
and regulate the activity of the integrin fil downstream effector 
molecule, integrin-linked kinase (ILK), which is known to play an 
essential role in GBM assembly and podocyte function [22]. In 
addition, ILK expression is stimulated in podocytes by renal injury 
and induces podocyte epithelial-mesenchymal transition [35]. 
Interestingly, the podocyte-specific knockout of ILK induces 
progressive focal segmental GS with increased mesangial matrix 
deposition and downregulation of Mphsl and Synpo with age, as 
observed in ICGN mice [12,36,37]. Activation of FAK, another 
effector molecule of the integrin pi signaling pathway, is also 
observed in the glomeruli of Tns2 mutant mice and conditional Ilk 
knockout mice [22,37]. This upregulation of FAK activity is 
considered to compensate for the ILK abnormality structurally, 



but not functionally [22,37]. It is notable that the Ilk knockout 
genetic background is derived from B6 {NPHS2 promoter-Cre 
transgenic mice) and 129 (ILK"™^™ mice) [36,37]. That is, the 
genetic backgrounds of B6 and 129 are resistant to Tns2- 
deficiency, but their resistance has no effect on //^-deficiency. 
Hence, it is assumed that the resistance effects are exhibited by 
functionally regulating the disorders in the integrin fil -ILK 
signaling pathway or structurally reinforcing the fragility of the 
actin cytoskeleton and slit diaphragm induced by aberrant ILK 
activity in Tn,s2-deficient podocytes. 

Under this assumption, we screened QTLs mapped on Chr 2 
for podocyte-related genes using XPodNet and identified 32 
candidate genes (Table 2). Among these, we were especially 
interested in Motchl, Dnml, Lmxlb, £eb2 and Wtl. NOTCH 1 
(18.91 cM) is a transmembrane receptor protein whose signaling 
pathway correlates with GS and renal function [38]. DNM1 
(22.09 cM), dynaminl, participates in endocytosis in the podocyte 
foot process and contributes to maintenance of the slit diaphragm 
[39]. WT1 (55.06 cM) and LMX1B (22.48 cM) are transcription 
factors that regulate podocyte-specific structural genes, Mphsl and 
Podxl, and JVphs2 respectively, and play a key role in the 
development and maintenance of the actin cytoskeleton and slit 
diaphragm in podocytes [40,41,42]. ZEB2 (27.31 cM), a tran- 
scriptional repressor of E-cadherin, is upregulated by the integrin 
P 1 -ILK signaling pathway downstream target SNAIL and induces 
the epithelial-mesenchymal transition [43]. 

In the same way, 1 4 candidate genes of the Chr 13 QTL are 
listed in Table 3. This locus partially overlaps the HIVAN 
susceptibility locus, HIVAN2, that was found in an (HIV-1 
transgenic FVB xB6) F 2 intercross [32]. An expression quantitative 
trait locus (eQTL) analysis has revealed that HIVAJV2 influences 
the transcript abundance of Mphs2, a membrane protein linking the 
slit diaphragm to the cytoskeleton [32]. However, HIVAN2 
includes none of transcription factors known to regulate Mphs2 
expression. In our cross, the locus was not linked to anemia, so that 
it is not considered to be essential for susceptibility to CKD 
induced by TnsT ph . 

Our genomic analysis of CKD in a cross between B6 and ICGN 
assumes that the Chr 2 loci critically contribute to resistance in 



PLOS ONE | www.plosone.org 



7 



June 2014 | Volume 9 | Issue 6 | e99602 



QTL for Resistance to Congenital Nephropathy in Tns2-Deficient Mice 



Tn.s2-deficient mice. The Chr 2 QTLs identified in crosses for 
proteinuria between resistant B6 and susceptible D2 or A/J or 
HIV-1 transgenic FVB are similar to ours [30,31,32]. This 
suggests that our QTL is not limited to resistance to the CKD 
induced by Tns2 nph , but is generally associated with resistance to 
CKD. Although this locus has not been noted in previous reports 
due to the absence of a homologous human QTL, we believe that 
this locus also plays an important role in resistance to CKD. 
Identification of resistant genes on this locus would help reveal the 
mechanism by which Tns2 deficiency leads to CKD and provide a 
novel insight into susceptibility to CKD, leading to the develop- 
ment of treatment strategies for CKD. 

Supporting Information 

Figure SI Representative examples of glomerular 
score. 

(TIF) 
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